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Abstract 

The physics regarding the existence of the critical end point (CEP) on the QCD phase boundary 
still remains unclear and its precise location is quite uncertain. In this paper we propose that the 
hard-core size of the baryons used in the description of the hot and dense hadron gas (HG) plays 
a decisive role in the existence of CEP. Here we construct a deconhning phase transition using 
Gibbs' equilibrium conditions after using a quasiparticle equation of state (EOS) for QCD plasma 
and excluded-volume EOS for the HG. We find that the first order transition results only when 
we assign a hard-core size to each baryon in the description of HG and the phase boundary thus 
obtained terminates at CEP beyond which a cross-over region occurs. The mean field approach 
for the HG lends support to this finding where unless we include an excluded-volume effect in 
the approach, CEP does not materialize on the QCD boundary. This investigation provides an 
intuitive reasoning regarding the origin of CEP and the cross-over transition on the QCD phase 
boundary. 
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I. INTRODUCTION 



Precise mapping of the QCD phase boundary existing between two distinct phases of hot, 
dense hadron gas (HG) and weakly interacting plasma of quarks and gluons (QGP) and 
the location of hypothesized critical end point (CEP) have emerged as interesting and chal- 
lenging problems before the experimental and theoretical heavy-ion physicists today [1-3]. 
The discovery of QCD critical end point is bound to clear the mist surrounding our under- 
standing of the conjectured QCD phase diagram and hence it would help us to ascertain the 
properties and signals of QGP [4-6] to some extent. The possible existence of CEP in the 
temperature (T) and baryon chemical potential (//s)plane of the QCD phase boundary was 
proposed a decade ago and it represents a second-order transition point where the first-order 
transition boundary terminates as T increases and fis decreases [3]. Its separation from the 
temperature axis (/i_B = 0) spans the region of a cross-over transition . Lattice QCD studies 
employing Monte Carlo simulation have failed at //b 7^ because the absence of a probabil- 
ity measure precludes direct computations and hence many mathematical approximations 
have been used to determine QCD phase diagram for nonvanishing values oi /ib- Some of 
these calculations predict the CEP to occur in the range jic/Tc ~ 1-0 — 1.4 [7-8]. However, 
certain calculations have also cast a shadow of doubt over the occurrence of a CEP in the 
phase diagram [9]. We still do not have any intuitive picture for understanding the cir- 
cumstances under which a cross-over transition can occur around /Ub = and which finally 
culminates into a CEP as //b increases. This is quahtatively supported by some lattice QCD 
findings [7-8]. In this paper, we take the help of a phenomenological model to emphasize 
the dominant role played by the finite-size baryons as constituents of a hot, dense HG in 
the existence of a cross-over as well as CEP on the QCD phase boundary. 
The search for a realistic equation of state (EOS) of HG is essential for the proper under- 
standing of the properties of the QGP. In a simple treatment of HG, all the baryons and 
mesons are treated as non-interacting point-like particles. However, such EOS of the HG 
has an undesirable feature that at very large and/or T, the hadronic phase reappears 
as a stable configuration in the Gibbs' construction of equilibrium phase transition between 
HG and QGP. Thus the pressure equahty condition Phg = Pqgp occurs twice for two values 
of /is for each T because of an exponential growth of hadrons and their resonances at higher 
^B [10]. Thus an anomalous feature of the reversal of phase transition from QGP to HG 



occurs which is indeed contrary to what we expect from the asymptotic freedom of QCD. In 
hterature, this feature was handled by modeUing the repulsive interactions existing between 
hadrons at large T and/or fiB of HG either in the mean field approach or in the excluded 
volume method. The attractive and repulsive interactions in the mean fields are incorpo- 
rated in HG by scalar a and vector w-exchange, respectively [11]. The Yukawa potential 
due to a;-exchange V{r) — [G"^ / i.T:r)exp{—m^r) generates a mean potential energy in HG 
as Ub — G^nB/rriij which vanishes when the net baryon density rie — >■ 0. This means that 
one can again generate a large number of hadronic resonances at large T where ub — and 
consequently Peg again becomes larger than Pqgp at very large T. In the recent past, we 
have attempted to cure this problem by adding a Vander-Waals repulsive interaction term 
UvDwi'f^i T) which depends on T and n (total number density of particles) and has its origin 
in the excluded volume correction [12, 13]. However, the main drawback of this model is that 
we cannot mathematically handle the EOS of HG if it incorporates many resonances in its 
description and a large uncertainty also results due to many unknown coupling parameters 
entering into the EOS. An alternate phenomenological description has mostly been used in 
the description of HG and the repulsive force arising due to hard-core volume of each baryon 
is geometrically incorporated as excluded volume correction in the pressure of HG [10, 14- 
16]. After such correction, we obtain a reduced pressure and hence Phg — Pqgp is satisfied 
for only one value of //^ at each T. In order to obtain Pqgp, we either use a simplistic bag 
model EOS or we take the help of more realistic quasiparticle model. Recently we pointed 
out that QCD phase boundary obtained from such exercise depicts a first order deconfining 
phase transition and the transition line indeed terminates at CEP beyond which a cross-over 
transition occurs when T further increases and fj,B decreases slightly [17-18]. 
In this paper, we plan to investigate the important role of baryon size as the origin in the 
existence of CEP and the resulting cross-over transition. We seek an answer to the ques- 
tion: why should a CEP occur on the conjectured QCD boundary and to what extent, the 
modeUing of the HG should be held responsible for the existence of CEP? 
The rest of the paper is organised as follows. In section 11, we discuss briefly the main 
features of the excluded volume models used in this paper. We also give outlines of our 
thermodynamically consistent formulation of the excluded volume model which was also 
used in the previous papers [17, 18]. Section HI deals with the mean field approach used 
for the description of HG where we incorporate an explicit term depicting the excluded vol- 



ume correction. In section IV, we give a quasiparticle model which is a thermodynamically 
consistent formulation of the EOS for QGP. Finally, section V gives our results, detailed 
comparisons and conclusions. 



II. FORMULATION OF EXCLUDED VOLUME MODEL FOR HG 

In the excluded volume models, repulsive interaction between two hadrons has been in- 
cluded by giving the hadrons a hard-core geometrical size and consequently reduces the 
hadronic degrees of freedom at large T and/or ^b- Consequently, the hadronic pressure 
is reduced and one can get a deconfinement phase transition from HG to QGP by using 
Gibbs' construction. However, some of these models arc not thermodynamically consistent 
because number density cannot be obtained directly from the partition function. There were 
many attempts to obtain a thermodynamically consistent formulation of the excluded vol- 
ume model [10, 14-16]. Recently we have proposed a thermodynamically consistent excluded 
volume model for hot and dense HG [17, 18]. Our approach has the following new features. 
Besides thermodynamical consistency, our model uses full quantum statistics so that the 
phase boundary in the entire (T, /i^) plane can be obtained without using any additional 
approximation. Recently, we have used our model for constructing a first order dcconfining 
phase boundary and the phase boundary is found to terminate at CEP [17, 18]. We have 
further determined the chemical freeze-out curve from our HG model and its proximity to 
CEP was pointed out [1] . In this model, we give importance to baryonic hard-core repulsion 
and we thus incorporate excluded-volume correction arising due to baryonic size only. We 
assume that mesons can overlap and fuse into one another and hence do not possess any 
hard-core repulsion. The grand canonical partition function for the HG, with full quantum 
statistics and after incorporating excluded volume correction in this model can be explicitly 
written as: 



where gi is the degeneracy factor of ith species of baryons, Ei is the energy of the particle 




(1) 



{Ei = ^Jk? + m^), is the eigenvolume of one baryon of ith species and Ylij ^j^j 



total occupied volume by the baryons and Nj represents total number of baryons of jth 
species. 

Now we can write Eq.(l) as: 



lnZr^V{l-Y,nfVf)k\, (2) 
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where Jj represents the integral: 

J _ 9i r k^dk 1 
'~G^Jo V^T^ [expif) + X,y 

and Xi — cxpIy) is the fugacity of the particle, n^^ is the number density of jth type of 
baryons after excluded volume correction and can be obtained from Eq.(2) as: 

This leads to a transcendental equation as 

nf = (1 - R)IX - /.A?|g + A,^(l - R)I[ (5) 

where /■ is the partial derivative of /j with respect to Aj and R — nf^V^ is the fractional 
occupied volume. We can write R in an operator equation as follows [10]: 

R^Ri + hR (6) 

where Ri = with R^ = + '^I'iV^^X^; n° is the density of pointlike baryons of 

ith species and the operator Cl has the form : 

"--T^E'WA.|- (7) 

i 

Using Neumann iteration method and retaining the series upto Cl'^ term, we get 

R^ Ri + nRi + n^Ri (8) 

Eq.(8) can be solved numerically for R. Finally, we get the total pressure [17, 18] of the 
hadron gas: 

p^H^a =T{1-R)J2 ^^^^ + E (9) 

i 3 

In (9), the first term in the right hand side represents the pressure due to all types of 
baryons where excluded volume correction is incorporated and the second term gives the 



total pressure due to all mesons in HG having a pointlike size. In this calculation, we 
have taken an equal volume — for each type of baryon with a hard-core radius 
r = 0.8 fm. We have taken all baryons and mesons and their resonances having masses upto 
2GeV/c^ in our calculation for HG pressure. Wc have also used the condition of strangeness 
neutrality by putting S'^nf — n|) = 0, where Si is the strangeness quantum number of the 
ith hadron, and is the strange (anti-strange) hadron density, respectively. We want 

to stress here that the form of this model used under Boltzmann approximation has been 
found to describe [21] the observed multiplicities and the ratios of the particles in heavy-ion 
collisions. In order to show the comparison of our results with the results obtained in a 
thermodynamically inconsistent approach of Cleymans and Suhonen [19] which has been 
described in detail in the ref. [10]. In this approach, the excluded baryon density of ith 
species can be written as [20] : 

III. MEAN FIELD MODEL FOR HG 

An alternate method for modelling the EOS for HG after incorporation of the repulsive 
interactions existing between hadrons of HG is the mean field approach. In this paper, we 
have used the mean field model of Tiwari et. al. [12], based on the work of Anchiskin 
and Suhonen [13], in which a Vander-Waals repulsive interaction term Uynwin, T) is added. 
This term has its origin in the excluded volume correction. We extend this model to describe 
the interactions in the HG and include the contributions of the baryons N, A, S and A- 
resonance in addition to the non-interacting mesons upto a cutoff mass of 2 GeV in HG. 
We have again treated the mesons as pointlike particles. The attractive interaction between 
the baryons is given by the scalar a-exchange while the exchange of vector w-meson gives 
the repulsive force. We have taken the value of coupling constant from one of the work of 
Suguhara and Toki [22] as follows 

= 3.178/m, — 3.871/m. (11) 

We use SU(6) quark model to obtain the relations [23-28] between various couplings as 

_2 _2 1 _ 

QujNN — -fi'tjAA, QujNN — gfl'wEE' '^9ojNN — dujEE- (12) 



Moreover, the couplings for A particle are assumed to be the same as those of nucleons [29] . 
The expression for the total pressure of the HG in this mean field model, with excluded 
volume correction, can be written as [12, 13] 

(13) 

In Eq. (13), the first term on the right hand side is the contribution from baryon with an 
effective mass M* and effective chemical potential /i*. Second term is the excess pressure 
because of the excluded- volume correction. Third term represents the baryon-density depen- 
dent mean-field pressure. Fourth term is the mean-field pressure due to cr-exchange and the 
last term on the right hand side is due to the contribution of pointlike mesons. Furthermore, 
the other terms involved in Eq. (13) are: 



i(i) 



, (Mf + ey/^ + UvDwin, T) ± UBjinB) T H \ , 
exp I I + 1 



(14) 



Here the upper (lower) sign refers to baryons(anti-baryons), respectively. The expression for 
Vander-Waal hard-core repulsion terms pvDW and Uynw are: 

PvDw{n,T)^nT-^^, (15) 
UvDwin, T) = - Tln{l - V^n), (16) 

1 — VqU 

with 

3 

and Vq is the hard core volume of each baryon (Vq = |7rrQ, tq = O.S fm). We have taken the 
same hard core volume for all type of baryons. Now the third term in the right hand side 
of Eq. (15) is represented as 

Pb = ^rnluj^, (18) 

where Uq is the time component of the w-exchange field and in the mean field approximation, 
it is given by 

ojQ^ [guiNNriBN + g^KKriBK + gojY: YPb e + gojssriBs + gojAAriBA] , (19) 

where 

/d^k 



nB = J2^Bj{fij,T). (21) 

j 

Here gj is the degeneracy factor. In order to calculate the net A number density, we use the 
relation [21] 

riBA ^n^-n^ = W{M)dM j d^k[U - /a]- (22) 

Here W{M) is the profile function which takes into account the finite width of A-resonance 
[30] . Furthermore, we get U bj in terms of the time component of vector field ou as 

UBjiriB) ^ Qojjji^o- (23) 

The vector interaction of o'-meson with all the baryons and A-resonance modifies their 
chemical potential as 

Pi* = pij - Ubj. (24) 
The attractive interaction of baryons with scalar field a modifies their masses as 

M* = M,- g^,,a, (25) 

and the pressure is 

Paia) = (26) 

In order to determine the effective mass of the different species, we have to determine the 
mean scalar field a. Using the thermodynamic consistency condition we can derive the 
following expression for the scalar field a as follows 



where 

[ d^k M* ^ 

^<rjQ) - 9j J ^2^^3 (^^2 + ^2)1/2 m (28) 

It is obvious to see that the solution for M* involves a set of coupled and self consistent 
equations and hence we have to solve the following set of seven coupled equations self 
consistently in order to get the values of the baryon densities and effective masses of different 



species. These equations are 



9A r r^. Ml 



+ sl^V (27r)3 (Mf + F)V2^^^ + ^a) (29) 
+ (27r)3(Mf2 + F)V2^^2: + /E)j 



1 /• d^k Mg 
3^" 7 (27r)3(M^2 + A;2)V2^-^= + -^^| 



and five equations for baryon density nBN,nBA,nBY,y''^BE,nBA as given by Eq. (20) and 
(22), respectively and one more equation represented by Eq. (17). Further, we get the 
effective masses of hyperons and A particles as follows 

^Xe = Ma,^-^(M^,-M^) 

= Ms - -{Mn - M^r) (30) 
3 

Ml = M^-{Mm-M%). 
For meson m, we have used the following ideal gas equation for its number density: 



In order to get the final result for Eq. (13), we impose the condition of strangeness neutrality 
to get the pressure of HG. It is obvious that we cannot include more baryons into the HG 
spectrum because the calculation becomes too much complicated to handle. 



IV. QUASIPARTICLE MODEL (QPM) 

We have used a thermodynamically consistent quasiparticle description as proposed by 
Bannur in order to study the EOS of QGP [31]. In this model, the system of interacting 
massless quarks and gluons can be effectively described as an ideal gas of "massive" non- 
interacting quasiparticles. The mass of these quasiparticles depends explicitly on T and 
implicitly on /ig via QCD running coupling constant. In this model, we start with the defini- 
tion of average energy and average number of particles and derive all the thermodynamical 



quantities from them in a consistent manner. The effective mass of the gluon changes with 
T and Hq as follows [32]: 



ml{T) = ^g\T)T^il+'-^], (32) 



6 



where Nc represents the number of colours. We have taken A^c = 3 in our calculation and: 



(33) 



Here Nf is the number of flavours of quarks and A*/ is the quark chemical potential belonging 
to the flavour f. Similarly the effective mass of the quarks involves the following relation 
[31]: 



rrii = mL + V2mgomth + 



m 



thi 



(34) 



■■"qO 

were rn^o is the rest mass of the quarks. In this calculation, we have used m^o = 8MeV 
for two light quarks {u,d), and m^o = SOMeV for strange quark. In the above Eq (34) rrith 
represents the thermal mass of the quarks and it can be written as 



8Nr 



7r2 



(35) 



Taking these values for the effective masses, energy density can be derived from the grand 
canonical partition function in a thermodynamically consistent manner and is given as [34]: 

-t4 oo 



1=1 



da 



ids 



(36) 



with ei{xil) = {xil)^Ki{xil) + 3{xiiyK2{xil), where Ki and K2 are the modified Bessel 
functions with Xi — ^ and index i runs for gluons, up-down quarks q, and strange quark s. 
Here di are the degeneracies associated with the internal degrees of freedom. Now, by using 
the thermodynamic relation e = — p, pressure of system at — can be obtained as: 



p{T,Hg = 0) 



Po 



+ 



dT 



e(T, Hg = 0) 

2^2 



(37) 



T To 

where po is the pressure at a reference temperature Tq. We have used po—^ at To=100 MeV 
in our calculation. We get the pressure for a system at finite ^g 



p{T,f,g)^p{T,0)+ / n, 

Jo 



d/i 



(38) 



TABLE I: Ratio of n^/n at the CEP using HG and quasiparticle model. 



HG Models coordinates of CEP r ub n UB/n 

(r,/XB) (/m) {fm-^) ifm-') 

Cleymans and Suhonen (166, 149) 0.8 0.163 0.86 0.190 

Our HG Model (166, 155) 0.8 0.104 0.54 0.192 

Mean field (163, 157) 0.8 0.0981 0.492 0.199 



where the expression for Ug is: 



oo ^ 

J2{-'^y~'pStnh{fijT)Ii{xil) 



(39) 



with Ii{xil) = {xilYK2{xil). Thus all the thermodynamical quantities can be obtained in a 
consistent way by using this model. We have used the following expression for the coupling 
constant [18] 

/ 



9\T) 

4:71 



(33 - 2Nf) In 
where = 115MeV and a = - 



At 



1 + a 



3 (153 - 19Nf 



In 21n^A/l + a 



(40) 



V. RESULTS AND DISCUSSION 



In Fig. 1, we have demonstrated the location of the CEP when we adapt QPM as the EOS 
of the QCD plasma and also use our excluded volume model for the description of HG. The 
open points 6,61,62 on the curve represent the variations in the coordinates of CEP when 
we change the hard-core radius of the baryons as r = 0.8/m, r = 0.6/m and r = 0.4/m, 
respectively Thus although CEP has its origin in the excluded volume effect, the drastic 
alterations in the hard-core radii do not yield much change in the location of CEP. Moreover, 
its proximity to the freezeout curve calculated in our HG model also remains unaltered. 

In Fig. 2, we have shown the variations in the phase boundary when we either change 
the EOS of the QCD plasma or we vary the EOS of HG. If we use simplistic bag model 
for the EOS of plasma and our excluded volume model for HG, we get the phase boundary 
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FIG. 1: QCD phase diagram mT — fiB plane. F is the freezeout curve obtained from our excluded 
volume model for HG. B is the first order deconfinement phase transition line using EOS in the 
QPM and EOS for HG in our excluded volume model. Open points b, bi and 62 are the locations 
of CEP for baryon's hard-core radius r = 0.8, r = 0.6 and r = 0.4/m, respectively. 

represented by the curve A and its terminal point a is the location of CEP. It lies closest to 
the freezeout curve. Similarly B, C and D represent the phase boundaries when we use QPM 
for the EOS of QGP, but the EOS of HG is respectively taken in the form of our excluded 
volume model, Cleymans and Suhonen excluded volume model, and mean-field model with 
excluded volume correction. Their terminal points b, c, and d represent the locations of 
CEP in these models, respectively. This exercise although results in a sizable variation in 
the coordinates of CEP, but the effect due to the details of the models is still found to be 
small. 

In Fig. 3, the details of the pressure equality condition Pqgp — Prg ^-nd around 
CEP are shown when we model QGP by QPM and HG by our excluded volume model. 
This clearly shows that this condition is very sensitive to a small variation of O.SMeF in the 




FIG. 2: QCD phase diagram mT — fiB plane. F is the freezeout curve obtained from our excluded 
volume model for HG. A is the phase boundary using bag model for the EOS of QGP [17]and 
our EOS for HG. The a gives the location of CEP. B is the boundary using EOS in QPM and 
our excluded volume model for HG and b gives CEP. C is the first order deconfinement phase 
transition line using QPM and the simple Cleymans and Suhonen model for HG and c is the end 
point of this curve. Similarly, D is the line obtained using QPM and mean field model for HG and 
d is the corresponding end point. 

temperature when we find that the condition fails to hold and the deconfining transition does 
not occur. It also defines the beginning of a cross-over region lying beyond the critical end 
point where the meson dominant HG pressure is always less than the QGP pressure. Hence 
the dual description with quarks, gluons together with tt and X-mesons as constituents 
applies well in this region. The main assumption in our model is that the pressure of 
HG is reduced even if a large number of hadrons is produced at large T and/or /ib- It 
means that at large ^b, the fractional occupied volume R in our model increases and hence 
mobility of the baryons is considerably reduced. However, in our consideration mesons do 
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FIG. 3: Variation of pressure with respect to //^ at different temperatures near CEP. Solid curve 
is the result for HG from our excluded volume model and dashed curve is the pressure for QGP 
obtained from QPM. Pressure equality condition is not fulfilled if T is increased by O.bMeV from 
T = 166MeV. 

not possess any such hard-core volume and hence they can fuse into one another when 
compressed. Thus we think of a possible parameter x which defines the ratio of baryon 
density ub to the total number density n{— ub + rim) in the HG at the CEP. Here rim is 
the meson number density. In table 1, we show the values of this ratio at CEP obtained 
in various prescriptions of HG we have adopted here. Surprisingly we find that the ratio 
X — 0.195 ± 0.005 which signifies that at CEP, the ratio Um '■ ub is almost fixed as 4 : 
1 and is independent of HG models used in the calculations. Beyond CEP, the meson 
density increases and thus yields more dominant contribution in the cross-over region. This 
demonstrates that the location of CEP on the phase boundary requires that all the baryons 
in the HG possess a hard core volume but mesons when being compressed, can fuse into each 
other. Thus ns/nm = 0.25 gives the location of CEP and this condition is independent of 
HG prescription used for the calculation. It still remains worth investigating problem why 
and how ns/rim = 0.25 yields the precise location of CEP. In the excluded volume model, 
this ratio defines a critical fractional occupied volume Rc- As /ib increases, the fractional 
occupied volume R by baryons increases and consequently the mobility of baryons in the 
hot and dense HG decreases fast and it finally results into a reduced pressure of HG so that 
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FIG. 4: Variation of normalized net quark density with respect to temperature at different fiq/T 
[35]. Different points are the lattice data from Ref. [36]. 

the Gibbs' conditions of equilibrium phase transition become satisfied. 

We must emphasize that we have used two different descriptions for QGP and HG, 
respectively. We find that the earher version of our model [21] for the HG used with the 
Boltzmann approximation describes well the ratios of various particle multiplicities and we 
hope that the present version with full quantum statistics will still improve the results of 
comparison with the experimental data. In order to have confidence in our EOS for QGP, 
we must test its predictions with the recent lattice results obtained at zero as well as finite 
baryon density. 

In Fig. 4 and 5, we have shown the variations of normalized net quark density Uq/T^ and 
normalized entropy density s/T^ with respect to T at different fig in the quasiparticle model 
(QPM). We find that our results yield a good fit to the lattice results. This comparison 
with the lattice calculation shows that QPM together with its parameters used here indeed 
gives a proper EOS for QGP even at finite In Fig. 6, we have attempted to show what 



p. /T=0.2 








■ ^ /T=0.4 

q 








|i /T=0.6 








1.1 /T=1.0 
q 

A 

A A 

A 


A 






A / 

/ 

/ 

A,' 








* 

' ▼ 

: : ' 








1 

1 

• ^ 










1 




— '^m^^ 






■ 


. \ .... \ ... . 




J 1 1 1 1 1 1 L 






FIG. 5: Variation of normalized entropy density with respect to temperature at /xb = 0.0 in QPM. 
Different points are the lattice data from Ref. [37] . 

happens to the change in the entropy density at CEP. We have calculated the difference 
^ = {s/T^)qgp — {s/T^)hg and demonstrated its variation with respect to the coordinates 
of the phase transition points lying at the boundary of Fig. 2. We find that ^ = 0.0 at 
the CEP and is minimum. Although we have not yet established that CEP obtained in our 
calculation is a second-order phase transition point. Most importantly it is the terminal 
point of the phase boundary. However, our results clearly indicate that it can either give an 
isentropic or a second order phase transition point. This is certainly an interesting finding. 
Although we have used two different models for the description of QGP and HG phases, the 
vanishing of net entropy density at the CEP outlines a continuity in these descriptions. 

The physical mechanism involved in this calculation is intuitively analogous to the per- 
colation model where also a first order phase transition results with 'jamming' of baryons 
and thus mobility of baryons is affected [38-39] . However, in the percolation model we do 
not have any comparison to what we should get in the QGP picture. Here we use a similar 




FIG. 6: Variation of {As/T^) = {s/T^)qgp — {s/T^)hg with respect to coordinates of various 
phase transition points on the {T,^b) phase boundary. We have used transition points from the 
curve B of Fig. 2. 

picture and we explicitly and separately consider both the phases, i.e., HG as well as QGP 
and hence it gives a clear understanding how a first-order deconfining phase transition can 
be constructed in nature and finally we reach an interesting finding that the baryonic size is 
crucially responsible for the existence of CEP on the phase boundary in such a construction. 
At low baryon density, overlapping mesons fuse into each other and form a large bag or 
cluster, whereas at high baryon density, hard-core repulsion among baryons, restricts the 
mobility of baryons. Consequently we consider two distinct limiting regimes of HG, one 
beyond CEP is meson-dominant regime and the other is baryon dominant region. 

A question generally arises : does our calculation offer any intuitive mechanism regarding 
the origin of the cross-over region? Beyond CEP, the cross-over region naturally appears in 
our model when T further increases and /is decreases. Cross-over is defined as a gradual 
change of the system from one phase to the other without a definite transition point. Lattice 



QCD has confirmed the existence of the cross-over region ai /ib — between HG and 
QGP. However, what happens in this region is still an open question so far as QCD is 
concerned. QCD involves two distinct vacua usually called as perturbative and physical 
one. Cross-over is thus realised by a gradual transition from one vacuum to the other. 
However, decomposition of hadrons to quarks and antiquarks one by one contradicts colour 
confinement because an isolated coloured object cannot exist in a physical vacuum [40]. 
In Nambu-Jona-Lasinio (NJL) as well as Polyakov extended Nambu-Jona-Lasinio (PNJL) 
models, the thermodynamic potential involves two degenerate minima at which two phases 
are in thermal, mechanical and chemical equilibria according to the Gibbs' criteria for the 
first order phase transition between the phases of broken and restored symmetry [41]. At the 
CEP temperature Tcep, the chiral transition changes to the second order. For T > Tqep, 
the thermodynamic potential has only one minimum and the transition is a smooth cross- 
over. However, the mechanism of cross-over is not understood very well because colour 
confinement does not strictly hold at the transition in these models. Our model falls in line 
with the ideas proposed recently [42-43] where it was shown that under circumstances, hot 
and dense HG consisting of extended hadrons could produce phase transition of the first or 
second order and also a smooth cross-over. We propose that although each baryon possesses 
a hard-core size, mesons are also extended particles but they lack a hard-core size. So they 
can overlap, fuse and interpenetrate. At CEP, mesons and baryons saturate the volume 
of the hot fireball. In meson dominated region (i.e., T > Tcep), niesons have a far larger 
density than that of baryons. When they start overlapping on each other, they fuse into one 
another and cluster formation starts where colour can fiow and only the cluster as a whole 
is colour-singlet. As the clusters merge together resulting into an infinitely sized cluster, 
analytic cross-over into a new phase occurs. Essentially we assume that each hadron is an 
extended bag of QGP and thus cluster formation arising due to fusion of mostly pions at 
T > Tcep, creates a smooth cross-over transition from one phase to the other. This picture 
appears more clear when we consider HG at iib — 0. However, we must emphasize that unlike 
other effective models, we use a hybrid model where EOS for HG and QGP are constructed 
independently and they reproduce separately the experimental particle-multiplicity data as 
well as lattice QCD results, respectively. It should be added here that many authors in the 
past have used two different equations of state for QGP and HG and obtained an explanation 
to an analytic and smooth cross-over and CEP in their models [44-45] . Our model presents 



a similar picture. Matching of the pressures at the CEP transition in the hybrid model used 
by us throws light on the mechanism of cross-over transition and as mentioned above, it is 
controlled by the presence of baryons in the system. But why does the ratio n^/nm = 0.25 
(a fixed value) at the CEP? How does the presence of baryon density affect the cluster 
formation? These questions need a thorough investigation before we make a clear picture. 

In conclusion, searching for the precise location of the critical end point (CEP) in the 
QCD phase diagram still poses a challenging problem. Although various calculations have 
predicted its existence but the quantitative predictions regarding its location wildly differ. 
Experiments face an uphill task in probing the CEP in QCD phase diagram because a 
clarity in theoretical prediction is missing. Moreover, many unstudied problems such as 
short lifetime and the reduced volume of the QGP formed at colliders also affect the location 
of CEP and its verification [46]. In these circumstances, our results arising due to baryon 
size, will be helpful in understanding the origin of CEP and determining its location on the 
phase diagram. 
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